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We present the first measurement of polarization and CP-violating asymmetries in a B^ decay 
into two light vector mesons, B^ — > (f>cj>, and an improved determination of its branching ratio 
using 295 decays reconstructed in a data sample corresponding to 2.9 fb~^ of integrated luminosity 
collected by the CDF experiment at the Fermilab Tevatron collider. The fraction of longitudinal 
polarization is determined to be /l = 0.348 ± 0.041(stat) ± 0.021(syst), and the branching ratio 
B{B^l-¥ 0<^) = [2.32 ± 0.18(stat) ± 0.82(syst)] x 10"^ Asymmetries of decay angle distributions 
sensitive to CP violation are measured to be Au = —0.007 ± 0.064(stat) ± 0.018(syst) and A^ — 
-0.120 ± 0.064(stat) ± 0.016(syst). 
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Several charmless 5° decays were observed at the 
Tevatron in Run II [3 |2], but a detailed investigation 
of decay properties and of CP violation in these decays 
is still lacking. The B^ — ^ ifxp process is mediated by a 
one-loop flavor-changing neutral current, the 6— >-s pen- 
guin, and belongs to the class of decays where the final 
state consists of a pair of light spin-1 mesons (V). Three 
independent amplitudes govern B — > VV decays, corre- 
sponding to the polarizations of the final-state vector 
mesons: longitudinal polarization, and transverse polar- 
ization with spins parallel or perpendicular to each other. 
The first two states are CP-even, while the last one is 
CP-odd. Polarization amplitudes can be measured an- 
alyzing angular distributions of final-state particles. In- 
terference between the CP-even and CP-odd amplitudes 
can generate asymmetries in angular distributions, the 
triple product (TP) asymmetries, which may signal un- 
expected CP violation due to physics beyond the stan- 
dard model (SM). 

The V-A structure of charged weak currents leads 
to the expectation of a dominant longitudinal polariza- 
tion [31 13] . Approximately equal longitudinal and trans- 
verse polarizations have been measured instead in 6— )■ s 
penguin-dominated B° and P"*" decay modes [5j . This is 
explained in the SM by including either non-factorizable 
penguin-annihilation effects [6] or final state interac- 
tions [7j. Recent theoretical predictions [21 U indicate 
a longitudinal fraction Jl in the 40-70 % range, when 
phenomenological parameters are adjusted to accommo- 
date present experimental data. Explanations involving 
new physics (NP) in the 6— >-s penguin process have also 
been proposed [5]. Additional experimental information 
in Bg penguin-dominated decays, such as B'^^ipip, may 
help distinguishing the various solutions [S], and can be 
used to derive upper limits for the mixing-induced CP 
asymmetries |10| . 

Triple product asymmetries are odd under time- 
reversal (T), and can be generated either by final-state 
interactions or CP violation. In flavor-untagged samples, 
where the initial B flavor is not identified, TP asymme- 
tries can be shown to signify genuine CP violation |llj . 
In this respect they are very sensitive to the presence of 
NP in the decay since they do not require a strong-phase 
difference between NP and SM amplitudes, as opposed 
to direct CP asymmetries [12 . The TP asymmetry is 
defined as Atp — r(TP>o)+r(TP<o) ' "^^ere F is the decay 
width for the given process. In B^^(p(j) decays two TP 
asymmetries can be studied, corresponding to the two in- 
terference terms between amplitudes with different CP. 
These asymmetries are predicted to vanish in the SM, 
and an observation of a non-zero asymmetry would be 
an unambigous sign of NP [12 . 

In this Letter we present the first measurement of 
polarization amplitudes and of TP asymmetries in the 
B^ — >■ (j)(p decay and an updated measurement of its 



branching ratio using P° -> J/ip(f> decays reconstructed 
in the same dataset as a normalization. Data from 
an integrated luminosity of 2.9 fb~^ of pp collisions at 
-/s = 1.96 TeV are analyzed. 

The components of the CDF II detector relevant for 
this analysis are briefly described below; a more complete 
description can be found elsewhere [T3]. We reconstruct 
charged-particle trajectories (tracks) in the pseudorapid- 
ity range jr^l < 1 [14j using a silicon microstrip vertex 
detector [15] and a central drift chamber [16|, both im- 
mersed in a 1.4 T solenoidal magnetic field. The detec- 
tion of muons in the pseudorapidity range \r]\ < 0.6 is 
provided by two sets of drift chambers located behind 
the calorimeters (CMU) and behind additional steel ab- 
sorbers (CMP), while the CMX detector covers the range 
0-6 ^ hi ^ 1-0 [TTj. A sample enriched with heavy-flavor 
particles is selected by the displaced-track trigger [IB] . 
based on the Silicon Vertex Trigger (SVT) [T5]. It pro- 
vides a precise measurement of the track impact param- 
eter {do), defined as the distance of closest approach to 
the beam axis in the transverse plane. Decays of heavy- 
flavor particles are identified by requiring two tracks with 
120 /^m < do ^ 1.0 mm and applying a requirement on 
the two-dimensional decay length, Pxy > 200 /xm [20] • 

We reconstruct B^ mesons by first forming (j)^K~^K~ 
and J/i/j — ^ IJ'^ ^^~ candidate decays from opposite-sign 
track pairs with mass within 15 and 100 MeV/c^ of the 
known [3T] (f) and J/'ip mass, respectively. At least one 
J/0 track is required to match a segment reconstructed 
in the muon detectors. We form 5° — )■ 00 (5° ~^ J/'4'4') 
candidates by fitting to a single vertex the (J/0 0) 
candidate pairs. In the P^— >■ J/00 case the fit constrains 
the mass of the two muons to the J/0 mass [H]. At 
least one pair of tracks in the P^ candidate must satisfy 
the trigger requirements. Combinatorial background and 
partially reconstructed decays are reduced by exploiting 
the long lifetime and relatively hard pt spectrum of P^ 
mesons. We follow closely the selection adopted in [T], us- 
ing the vertex fit x'^, the Lxy, the reconstructed B° and 
meson impact parameters, and the minimum kaon trans- 
verse momentum as discriminating variables. The se- 
lection requirements are set by maximizing the quantity 
S/\/S + B, where the accepted number of signal events 
S is derived from a Monte Carlo (MC) simulation [52] 
of the CDF II detector and trigger, while the number of 
background events B is modeled using data in mass side- 
band regions: (5.02,5.22) and (5.52,5.72) GeV/c^. The 
resulting mass distributions are shown in Fig. IT] 

A binned maximum likelihood (ML) fit to the ttib 
distribution is performed to determine the P^ yield for 
both decay modes. The signal is parameterized by two 
Gaussian functions with the same mean value, but dif- 
ferent widths. The ratios between the two widths and 
between the integrals of the two components are fixed 
based on MC simulations. The combinatorial back- 
ground has a smooth mass distribution near the signal 
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FIG. 1: The invariant mass of the four kaons (left) and of 
the J/ij) and two kaons (right) for B^ — > (f><f) and B^l -^ J/'<P4> 
candidates, overlayed with fit projections and separate signal 
and background components. The narrower signal peak for 
the -Bs — >■ J/ip4> is due to the J/ip mass constraint applied in 
the reconstruction. 



and is modeled with an exponential function. A reflec- 
tion from B" -> </>X* (892)0 (^o ^ J/^K*(892)0) with 
misassigned kaon mass to final state pions contaminates 
the B^ — >■ 4)(j) (i?° — >■ J/ipcf)) signal region. Parameteri- 
zations and efficiencies determined from simulation are 
used for these backgrounds. Their normalizations are 
derived from the known |21| branching ratios, fragmen- 
tation fraction ratio fs/ fd, and the ratio of the detec- 
tion efficiencies relative to signal ones. We estimate 
(4.19 ± 0.93)% and (2.7 ± 1.0)% reflection background 



under the i?° — >• J/00 and B" 



-^ 



signals, respec- 



tively. Free parameters of the fit are the signal frac- 
tion, the _B^ mass M, and width a, together with the 
exponential slope &o defining the combinatorial back- 
ground mass shape. We estimate the total number of 
signal decays as iV^^ = 295 ± 20(stat) ± 12(syst) and 
Njp^ — 1766 ± 48(stat) ± 41(syst), where the system- 
atic uncertainty is estimated by varying signal and back- 
ground models. 

The _B^— >00 decay rate is derived from the relation 

g(Bg^00) ^ iV^ejJAA^w^e^ ^ 
B {B^,^J/U) N^^ B{cj,^ KK) e^^ ''^^ ' 

where e^^/e^^ is the acceptance times efficiency ratio for 
the two decays and ejl j, is the efficiency for identifying 
at least one of the two muons. The efficiency ratio is 
determined using a MC simulation of the CDF II detec- 
tor and trigger, whose reliability in determining relative 
trigger and reconstruction efficiencies has been verified 
for several different decay modes also using data-driven 
approaches |23j. We estimate e^^^/e^^ = 0.939 ± 0.099, 
where the uncertainty includes systematic effects from 
polarization uncertainties in the two decay modes (9%), 
from the different trigger efficiencies for kaons and muons 
(4%), and from the S" pr spectra (1%). We use in- 
clusive J/0 data to derive the single-muon identifica- 
tion efficiency as a function of muon px- It is deter- 
mined separately in two pseudorapidity regions corre- 



sponding, respectively, to the CMU/CMP and CMX de- 
tectors, and is described by a turn-on function that de- 
pends on a plateau, a slope, and a threshold parame- 
ter. We use simulated i?° —> J/00 decays to calculate 
ejt , treating the efficiencies for the two muons as uncor- 
related: e^^;^ = (86.95 ± 0.44(stat) ± 0.75(syst))%. The 
systematic uncertainty includes the uncertainty on the 
background subtraction and effects of residual correla- 
tion between the two muon efficiencies. 

We measure S(B^ -^ 4>4>)/B{B° -^ Jj-^c^) = [1.78 ± 
0.14(stat) ± 0.20(syst)] x lO'^ and derive S(B°-^00) = 
[2.32±0.18(stat)±0.26(syst)±0.78(br)] X 10~^ using the 
known [21, B{Bfl^ Jl-\l}(l))^ which contributes the domi- 
nant uncertainty, labeled (br). This result is in agree- 
ment and supersedes our previous measurement [!• with 
a substantial reduction of its statistical uncertainty; it is 
also consistent with recent theoretical calculations [31 S] . 

We describe the angular distribution of the 5° -^ 
00 decay products using the helicity variables ui — 
(cos^i,cosi?2,^), where '&i is the angle between the di- 
rection of the K^ from each and the direction op- 
posite the i3° in the vector meson rest frame, and <? 
is the angle between the two resonance decay planes 
in the 5° rest frame. The three independent complex 
amplitudes are Aq for the longitudinal polarization and 
A^ (^_l) for transverse polarization with spins paral- 
lel (perpendicular) to each other. They are related by 
|Aop -H |A||p + |yl_Lp = 1. The differential decay rate is 

expressed as d'^T / (dtdoj) ex J2i=i ^ii^) fii^) ^ where the 
functions Ki{t) encode the i?" time evolution including 
mixing and depend on the polarization amplitudes, and 
the fi{ijj) are functions of the helicity angles only |12) . 
To extract the polarization amplitudes we measure the 
time-integrated angular distribution assuming no direct 
CP violation and a negligible weak phase difference be- 
tween _B^ mixing and Bl^xfxj) decay as predicted in the 
SM. The time-integrated differential decay rate depends 
on the polarization amplitudes at t = and on the light 
and heavy B'^ mass-eigenstate lifetimes, tl and th, as 
follows: 

d^V 

-^ (X rL(|Aop/iM + |A||P/2(^) 

+ |Ao||A|||cosJ||/5(^)) +rH|AiP/3(^), (1) 
where (5|[ — arg(AQA||) and 

fi{uj) =4 cos^ d 1 cos^ i?2 , 
/2(w) =sin2i?isin2t?2(l + cos2<?), 
/3(tj) =sin2i?isin2t?2(l-cos2<?), 
f^iuJ) = y2sin2i9isin2i92COS<?. 

Two triple products are present in _B — >■ VV decays: 
TP2 = 3(Aj|^_l), and TPi = '^{AIAa_). These factors 
appear, respectively, in the decay rate terms 1^4 (t) and 



KQ{t) multiplied by the functions 

/4(a;) = -2 sin^ i?i sin^ iJj sin 2$, 
/g(cj) = -y2sin2?9i sin2^2sin^. 

In flavor-untagged samples tlie TP terms, that vanish 
in the absence of NP, are proportional to the so-called 
true triple products, and provide two CP-violating ob- 
servables, ^^p ^^^ A^p [H]. We access A^p through 
the observable u — sin 2^. We measure the u asymme- 
try, Au, by integrating over cos'i?i,2 the untagged decay 
rate and counting events with u > (iV+) and u < 
{N~). Similarly, A\p is accessed through an asymme- 
try in sin^. We define the observable w as w = siwP 
{v = — sin^) if cosz9i cos??2 > (cosi^i cosi?2 < 0) and 
measure its asymmetry Ay by counting events with v>0 
(N.^) and v<0 (N^). The asymmetries are defined as 

A — "(") "(^) _ A/- „ 



N' 



u\v) u{v) 



[3(A]r.^^A 



where the two normalization factors are J\fu = — 2/7r and 
Afy = —\f2l'K. Both Au and Ay are proportional to 
CP-violating TP asymmetries, and are also sensitive to 
mixing-induced TP when considering the decay-width 
difference of the _B° system. 

We perform an unbinned ML fit to the reconstructed 
mass of the _B" candidates and the helicity angles in or- 
der to measure the polarization amplitudes. The con- 
tribution of each candidate to the likelihood is Li — 
fs'PsimBi,uJi\£.s) + (1 - fs)'Pb{mBi,uJi\^b), where /« is 
the signal fraction and Vj are the probability density 
functions (PDFs) for the B^ — ?> (/>(/) signal (j = s) and 
background (j = b) components, which depend on the fit 
parameters ^s and ^f,, respectively. The effects of neglect- 
ing the reflection background are included in the system- 
atic uncertainties. Both the signal and the background 
PDFs are the products of a mass component, described 
earlier, and an angular one. The signal angular com- 
ponent is given by Eq. [T] multiplied by an acceptance 
factor. The acceptance is computed in bins of the he- 
licity angles from simulated iJ^ — ^ </>(/) decays averaged 
over all possible spin states of the decay products and 
passed through detector simulation, full reconstruction, 
and analysis cuts. We use an empirical parameterization 
derived from the observed angular distributions in the 
mass sidebands to model the background angular PDF: 
the product of a flat distribution for the ^ angle and a 
parabolic function for the other two, whose single pa- 
rameter bi is a fit parameter. We fix tl and th to the 
world average values j21j . There are eight free param- 
eters in the fit: fs, ^s = (-^iCrJ^oPi l^||P: cos Jy) and 
Cb — {bo,bi). The fit has been extensively tested using 
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FIG. 2: Angular distribution for B^^<f>4> events with the fit 
projection, signal, and background component superimposed. 



simulated samples with a variety of input parameters and 
shows unbiased estimates of parameters and their uncer- 
tainties. We also perform the polarization measurement 
using the sample of «1700 B^ — > J/00 candidates de- 
scribed earlier. We find |^oP = 0.534 ± 0.0I9(stat) and 
|A||p — 0.220 ± 0.025(stat), in good agreement with cur- 
rent measurements |24) . The results of the polarization 
analysis for the B^^xjxj) sample are summarized in Ta- 
ble [ij In Fig. [2] we show the fit projections onto the helic- 
ity angles. The dominant correlation of the fit parame- 
ters is between |^oP and |^||P (-0.447), the others being 
much smaller. Several sources of systematic uncertainty 



TABLE L Summary of the 5° — >■ 4'4> measurements. The first 
uncertainty quoted is statistical and the second is systematic. 



Observable 



Result 



B 

cos(5|| 

Ay 



[2.32 ±0.18 ±0.82] x 10"^ 
0.348 ±0.041 ±0.021 
0.287 ±0.043 ±0.011 
0.365 ±0.044 ±0.027 

-0.9lj:[i:}^±0.09 
-0.007 ±0.064 ±0.018 
-0.120 ±0.064 ±0.016 



have been studied. We account for the neglected physics 
backgrounds considering the i3°— )-0iir*(892)° decay and 
two other possible contaminations: 5°— >-0/o(98O), with 
fQ-^K+K-,an.dB°^(j)K+K- (non-resonant). The lat- 
ter two contributions are normalized to the signal yield 
in analogy with similar B^ — ^ 4>X decays. We assume 
up to 4.6% contamination from B^ — >■ (/)/o and 0.9% of 
B° -^ cj)K+K-, and determine a 1.5%(0.4%) shift in 
the central value for |ylop(|A||p) using simulated exper- 
iments. Biases introduced by the time integration are 
examined with MC simulation: they are created by the 
dependence of the angular acceptance on AFs and by 
a non-uniform acceptance in the 5° proper decay time 
introduced by the displaced-track trigger. The assigned 
systematic uncertainty (1%) is the full shift expected in 
the central value, assuming a value for AF^ equal to the 
world average plus one standard deviation [ST. We also 



consider the propagation of ti^(h) uncertainties to the po- 
larization amphtudes (1%). Other sources of minor sys- 
tematic uncertainties are the modeling of the combina- 
torial background (0.4%) and of the angular acceptance 
(0.5%). The impact of CP-violating effects on the mea- 
sured amplitudes is negligible. 

The asymmetries Ai {i = u, v) are evaluated through 
an unbinned ML fit to m^ only, using the joint like- 
lihood for the N^ and N[' events with positive and 
negative u{v). The same tob PDF parameterization 
discussed above is used for samples with both u{v) 
signs. We multiply the total likelihood by the binomial 
f{N^,N~\p), where the probability p of obtaining N^ 
and N~ events depends on the overall signal fraction fs, 
the signal asymmetry Ai, and the background asymme- 
try A* : p = i [1 + AJs + {l~ fs)Al] . Mass and width for 
the i?g signal, as well as signal fraction, are consistent 
with those obtained in the polarization analysis, while 
background asymmetries are consistent with zero. The 
measured B^^KJxf) asymmetries are reported in Table H] 
The systematic uncertainty is evaluated using an alter- 
nate background parameterization as in the polarization 
analysis and by conservatively assigning maximal asym- 
metry to the neglected physics background peaking in the 
signal region. Using a large sample of simulated events, 
we check that the detector acceptance and resolution in- 
troduce a bias in the asymmetries smaller than 0.2%. 

In summary, we measure for the first time the polar- 
ization amplitudes and the triple product asymmetries 
in the B^^ffxp decay. We find a significantly suppressed 
longitudinal fraction /l = |Aop = 0.348 ± 0.041(stat) ± 
0.021 (syst), smaller than in other h-^s penguin B^VV 
decays [5j. This result agrees well with predictions [5] 
based on QCD factorization, but only marginally with 
perturbative QCD ones [3], and hints at a large penguin 
annihilation contribution [9]. The two measured asym- 
metries are statistically consistent with the no CP viola- 
tion hypothesis, although Ay is 1.8 cr different from zero. 
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